
%  fct_s =@ (x) exp(-(1./((x-1).^2) + 1./((1-x).^2)));
%  
%  fct_nu=@ (x) zeros(size(x))+(x>=1) + (x<1).*(x>0).* ( fct_s(x-1) ./ ( fct_s(x-1) + fct_s(x)));

%fct_nu=@ (x) x;

%fct_nu=@ (x) 3*x.^2 -2*x.^3;

fct_nu=@ (x) 5*x.^3 - 5*x.^4 + x.^5;

fct_W=@ (x) zeros(size(x)) + (x<=5/6).* (x>= 2/3).* cos(pi/2 .*fct_nu(5-6.*x))...
                                  + (x<=4/3).* (x>5/6)...
                                  + (x<=5/3).* (x>4/3).* cos(pi/2.*fct_nu(3.*x-4));
                              
fct_V=@ (x) zeros(size(x))+ (abs(x)<=1/3) + (abs(x)>1/3).*(abs(x)<=2/3).*cos(pi/2.*fct_nu(3.*abs(x)-1));

x=-2:.01:2;
y=0:.01:2;

figure;
my_ax=axes('fontsize',13);
plot(x,fct_V(x),'LineWidth',2);
set(my_ax,'XTick',[-0.66;-0.5;-0.33;0;0.33;.5;0.66],'XTickLabel',{'-2/3','-1/2','-1/3',0,'1/3','1/2','2/3'});
set(my_ax,'YTick',[0;.5;1]);
set(my_ax,'XGrid','on','YGrid','on');
axis([-.8 .8 0 1.05]);
xlabel('\omega','fontsize',20);
ylabel('V(\omega)','fontsize',20);


% figure;
% my_ax=axes('fontsize',13);
% plot(y,fct_W(y),'LineWidth',2);
% set(my_ax,'XTick',[.6667;.8333;1.3333;1.6667],'XTickLabel',{'2/3','5/6','4/3','5/3'});
% set(my_ax,'YTick',[0;.5;1]);
% set(my_ax,'XGrid','on','YGrid','on');
% axis([.5 1.833 0 1.05]);
% xlabel('r','fontsize',20);
% ylabel('W(r)','fontsize',20);

% figure;
% my_ax=axes('fontsize',13);
% plot(y,fct_W(y),'LineWidth',2);
% set(my_ax,'XTick',[.5;.6667;.8333;1;1.3333;1.6667;2],'XTickLabel',{'1/2','2/3','5/6','1','4/3','5/3','2'});
% set(my_ax,'YTick',[0;.5;1]);
% set(my_ax,'XGrid','on','YGrid','on');
% axis([.5 2 0 1.05]);
% xlabel('r','fontsize',20);
% ylabel('W(r)','fontsize',20);